DOUBLING PROPERTIES OF CALORIC FUNCTIONS 
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"Science is the great antidote to the poison of enthusiasm and superstition" (Adam Smith, 

Wealth of Nations, 5, III, 3) 



Abstract. We obtain quantitative estimates of unique continuation for solu- 
tions to parabolic equations: doubling properties and two-sphere one-cylinder 
inequalities. 



1. Introduction 

This work is concerned with quantitative estimates of unique continuation for 
parabolic equations. In the context of elliptic equations the quantitative estimates 
of unique continuation, which have been derived from the nowadays standard meth- 
ods to prove qualitative results of unique continuation ( Carleman or frequency func- 
tion methods) are the so called doubling properties and optimal three spheres in- 
equalities. In particular, it is by now well known that when 

n 

i X WW- ) 

is an elliptic operator in R™, there is a constant N depending on n, M, the parame- 
ters of ellipticity and the Lipschitz norm of the matrix of coefficients of the operator 
E such that, solutions of the inequality, \Eu\ < M (\u\ + \Vu\) on B4, verify the 
doubling property [11] 

u 2 dx < (NQ) N [ u 2 dx , 

B 2r (z) JB r {z) 

when \z\ < 1, < r < 1/2 and 

6 = / u 2 dx/ u 2 dx 

J £?4 J B\ 

and the optimal three sphere inequality ([13j. [14] ) 

1 N log (1/ r) 

ft \ 1 + Nlogd/r) / f \ l + Nlog(l/r) 

u 2 dx < { u 2 dx] \ N I u dX 



' B\ \J B r / \ J B^ 

when < r < 1. 

The three sphere inequality is called optimal for the following reason. If one 
seeks to find the largest possible function : (0, 1) — > (0, 1) such that the three 
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sphere inequality 





) 



8(r) 




) 



l-e(r) 



(1.1) 



u 2 dx < 



u 2 dX 



holds for some positive constant N, for all harmonic functions u in and < r < 1, 
one finds from the identity 



which holds when u is a homogeneous harmonic polynomial of degree k > 1, and 
(fj~Tj) that 



Thus, (9(r) cannot be larger than log4/log (4/r). On the other hand, the recent 
applications to stability issues in Inverse Problems of quantitative estimates of 
unique continuation ([2], [3]) show it is useful to know that "harmonic" functions 
are as small as they can be at intermediate scales when they are known to be small 
at smaller scales. 

The first quantitative result of strong unique continuation for parabolic equa- 
tions, a two-sphere one-cylinder inequality [13[ Theorem 9'] (see the body of the 
paper for the relevant definitions), was derived in 1974 from the corresponding first 
qualitative results of strong unique continuation for second order parabolic equa- 
tions in the literature [13]. In [13], E.M. Landis and O.A. Oleinik used a reduction 
argument to the already established elliptic results of unique continuation, which 
relays on the representation formula of solutions to parabolic equations with time 
independent coefficients in terms of the eigenfunctions of the corresponding elliptic 
operator, and thus their results did not apply to parabolic equations with time 
dependent coefficients. After 1988 when more results of unique continuation for 
elliptic equations were available in the literature, the qualitative results in [13j were 
reproved with less regularity assumptions in [14] . The quantitative results in [13] 
were also reprove with less regularity assumptions in [5j. The authors in [14] and 
[5] used the same reduction argument as in [13] and as a consequence, their results 
only apply to parabolic operators with time independent coefficients. 

In the same way as it has not been a trivial task to derive quantitative estimates 
of unique continuation for elliptic operators from the corresponding qualitative 
results (See [11] and [14]), the same happens in the parabolic case. 

In [16] , [6] , [7] , [8] , [9] , [10] and [1] the authors have obtained with standard argu- 
ments of unique continuation, frequency functions or Carleman methods, qualitative 
results of strong unique continuation for parabolic equations with time dependent 
coefficients. In particular, between [|] and [TO] the following qualitative-quantitative 
result is proved in relation with backward parabolic operators 




1/4 < (r/4) e for all k > 1 . 



n 

P = E 9 i (g ij (X)d j )+d t , 
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which for some positive constants A and M satisfy the conditions 

n 

(1-2) A|C| 2 < V.9 y '(X)^<A- 1 |^| 2 , 



,1/2 



(1.3) ]T \g^(X)-9 i HY)\<M(\x-y\ 2 + \t-s\y 

when X = (x, i), F = (y, s) are in R" +1 and £ is in M™. 

Theorem 1. Let P be a backward second order parabolic operator as above. Then, 
there is a constant N = N(X, M, n) such that, if u satisfies 

\Pu\ < M{\u\ + |Vu|) 

in Q2 and u(x, 0) has a zero of infinite order at x = 0, the following holds in Q\ 

\u(x,t)\<Ne-^\\u\\ Loa{Q2) . 

Here and in the sequel B r (z) — { x £ l n : \x — z\ < r }, B r — B r (Q), 
Q r {z, t) = B r (z) x [r, t + r 2 ] and Q r — Q r (0, 0). 

Theorem [T] gives a qualitative-quantitative result of strong unique continuation 
but is missing the proper and natural quantitative estimates of strong unique con- 
tinuation for second order parabolic equations: doubling properties within char- 
acteristic hyperplanes and two-sphere one-cylinder inequalities. Here, we prove 
these quantitative estimates with new arguments, which are based on the frequency 
functions or Carleman inequalities appearing in [16j . [3] and [10] and extend their 
applicability to solutions of parabolic equations with time dependent coefficients. 
In particular, the following theorem is proved: 

Theorem 2. Let P be a backward second order parabolic operator verifying (|1.2p 
and (|1.3|) and u satisfy \Pu\ < M (\u\ + |Vu|) in Q4. Then, there is a constant 
N = N(X, M, n) such that the following properties hold when < r < 1/2 

(1) J B2r U 2 (x,0)dx < (NQ) N J Br U 2 (x,0)dx. 

W log (!/■■) 



(2) f Bi u\x, 0) dx < (j Br u\x, 0) dx) 1+WtosWrJ (N ! Qi u 2 dx) 



l + JVlog(l/r) 



(3) Ifu\ t= o is not identically zero, then \u\l =0 is a Muckenhoupt weight in Bi/ 2 - 
Here, Q = Jq^ u 2 dX/ J Bi u 2 {x, 0) dx and dX — dxdt. 

Recall that a nonzero and locally integrable function w in B\ is called a Muck- 
enhoupt weight when there are numbers p £ (1, +00) and C > such that 




< C-h w dx and / w dx < C I 

Jb,.(z) JB 2rM JE 



w dx 

'BrW 



for all balls B r (z) such that B± r (z) C B\. Property 2 generalizes the analogue 
result for elliptic operators in [11] . 

In regard to space-time like doubling properties we have the following result. 

Theorem 3. Assume that P and u satisfy the conditions in Theorem® Then, there 
is N = N{X,M,n) such that the following holds when < r < l/y/N log (NO) 

(1) J Qr U 2 dX < e Nlog(N0)lo S (Nlog(N0)) r 2 f g ^ X)0 ) dx. 

(2) J Q [ r U 2 dX < e Nlo S (N&)lo S (Nlo S (N&)) J Q J u 2 dX _ 
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We do not know if Theorem [3] is optimal. In fact, if u is a backward caloric 
polynomial of degree k > 1, (i.e. u(Xx, X 2 t) = X k u{x,t) when (x,t) is in R n+1 , 
A > 0) the following happens when r > 



r 2 J B u 2 (x,0)dx J B u 2 (x, 0) dx 

Theorems [5] and [3] also hold when u verifies \Pu\ < M (\u\ + |Vu|) in Q 4 n D x 
[0, +oo), D = { (x',x n ) : x n > tp(x') }, where ip : M n_1 — > K is a C 1 ' 1 function 
verifying ip(0) — and when either the Dirichlet or Neumann (conormal derivative) 
data of u vanishes identically on dD x [0, +oo) n Q4. In fact, when the matrix of 
coefficients of P is time independent and u has zero Dirichlet condition on the 
lateral boundary the results hold when <p — ipi + <p2, where ipi is a convex function 
and (fi2 is C 1,a for some a in (0,1). Of course, one must replace B r by B r n D 
and so on in the corresponding statement of the boundary version of the Theorems. 
We do no prove these results here but they follow combining the arguments in this 
work with others in [9] . 

In section [2] we prove Theorems [2] and [3] when P = A + dt with the frequency 
function method. With some extra work it is possible to prove Theorems [2] and [3] 
with the same method and when the main coefficients of P are time independent. 
In section [21 the more general case is proved using a more detailed version of the 
Carleman inequalities appearing in [9] and [10] . 

2. Constant coefficients 

The proof of Theorems [2] and [3] when P = A + dt is a consequence of the following 
four Lemmas. 

Lemma 1. Assume that u satisfies \Au-\-dtu\ < M (|u| + |Vu|) in Qi. Then, 
there is N — N(n, M) such that the inequalities 



AHog(iVe) > 1 and N I u 2 (x,t)dx> / u 2 (x,0)dx 

JB 2 JBi 

hold when t < 1/ (iVlog(JV8)) and 9 = J Qi u 2 dX/ J Bi u 2 (x, 0) dx. 

Proof. Let / = u<p, where <p S C£°(B2), < ip < 1, tp = 1 in -83/2- The standard 
estimate [H] 

(2- 1 ) lklU»(Q 3 ) + l|V«|U«. (Qa ) < N\\u\\ L 2 (Qi) 

for solutions of the inequality 

(2.2) \Au + d t u\ < M (\u\ + |Vu|) in Q 4 
gives 

(2.3) |A/ + d t f\ < N [|/| + |V/| + \\u\\ LHQi)X B 2 \B 3/2 

and shows that the first claim holds. Setting H(t) = J f 2 (x, t)G(x — y, t) dx, where 
G(x, t) = t-™/2 e H*l 2 /4t and y e Biwe have 

(2.4) H(t) = 2 / /(A/ + d t f)G(x - y, t) dx + 2 / | Wf\ 2 G(x - y, t) dx 
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and from (|2.3j) . (12. 4|) and the Cauchy-Schwarz's inequality 

Integration of this inequality in (0, t) for t € (0, 16) gives 

N J f(x,t)G(x-y,t)dx>u 2 (y,0)-Ne- 1 / m \\u\\l 2{Qi) . 
Integrating this inequality over B\ and recalling that J G(x — y,t) dy = 1 we get 

Nj u 2 {x,t)dx> j u 2 (x,0)dx - Ne- 1/Nt \\u\\\ 2(Qi) , 
and the Lemma follows from this inequality and (|2.1|) . □ 
Lemma 2. Given a > and f £ W 2 ' 00 {Rl +1 ) set G a = (t + a)~ n/2 e -M 2 /4(t+a)^ 

H a (t) = f f 2 G a dx , D a (t) = f I V/| 2 G Q dx and N a (t) = 2> ' ~ " )LK ' ' ' 

JR™ JR™ 



H a {t) 



Then, 



N a (t) > J (A/ + d t f) 2 G a dx 



Proof. The identities d t G a - AG a = 0, VG a = -j^G a , A = div(V ) and 
integration by parts imply the following identities 

(2.5) H a (t) = 2 J f(Af + d t f)G a dx + 2D a (t) , 

H a (t) = 2 J f(dtf + ^a) -V/-i(A/ + a t /)) G a dx + J f(Af + d t f)G a dx, 

D a (t) = J f (d t f + ■ V/ - § (A/ + dtfj) G a dx-\ J f(Af + d t f) G a dx 

and 

(2.6) H a (t)D a (t) = 2(J f(d t f+^-\7f-±(Af + d t f))G a dx^ 

-\(Jf{Af + d t f)G a dx^ . 

The Rellich-Necas identity with vector field VG a 
div(VG a |V/| 2 ) - 2div((V/ ■ VG a )V/) 

= AG a |V/| 2 - 2D 2 G a Vf ■ V/ - 2V/ ■ VG a Af 
and integration by parts give 
(2.7) 

J AG a \Vf\ 2 dx = 2 J D 2 G a Vf- Vfdx + 2 J Vf-VG a Afdx 

= 2 / (^-V/) 2 G a dx-2 J ^.WfAfG a dx-D a (t)/(t + a) . 

Again, the fact that G a is a caloric function, integration by parts, (|2.7|) and the 
completion of the square of dtf + 2 (t+a) f ~ § (A/ + c?t/) yields the formula 



6 



L. ESCAURIAZA, F.J. FERNANDEZ, AND S. VESSELLA 



(2.8) D a (t) = 2 J (a t /+^-V/-i(A/ + 9 t /)) 2 G Q da: 

- I / (A/ + ft/) 2 G a dx - D a (t)/(t + a) . 
Then, from (|2l) ]) .([2~%] ) and the quotient rule 

(2-9) N a (t) = ^±^U (d t f + ^-Vf-±(Af + d t f)YG a dxH a (t) 

- (/ / (ft/ + ^ ■ V/ - | (A/ + ft/)) G a ds) 

+ \ (^J f (A/ + ft/) G a dx) - \ J (A/ + ft/) 2 G a dxff a (i) j 

and Lemma [5] follows from (|2.9p , the Cauchy-Schwarz inequality and the positive- 
ness of the third term on the right hand side of (|2.9[) . □ 

Lemma 3. The inequality 

J \^ e -M 2 /ia dx < 2a J |V/i| 2 e-l K l 2 / 4a dx + f J h 2 e-^ 2 ' ia dx 
holds for all h £ G °°(M n ) and a > 0. 

Proof. The inequality follows setting v = he~\ x \ / 8a and from the identity 
2a J \Vh\ 2 e~^/ 4a dx+^ J /i 2 e -M 2 / 4 ° dx - J JfJ^e-N 2 ^ dx 



= 2a J \Vv\ z dx . 

□ 

Lemma 4. Assume that N and 6 verify Nlog(NO) > 1, h e Cg°(R n ) and that 
the inequality 

2a / 1 * + f / d. < » log (WO) / ^ * 

ZioWs w/ien a < 12 at i g (Are) • Then, 

h 2 dx < (NQ) N / ft 2 da; w/iera < r < 1/2 . 
Proof. The inequality satisfied by h and Lemma [3] show that 

y J|£/i2 e -l«.| a /4« dx <N\og{NQ) J h 2 e~\ x \ 2 ' ia dx 
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when a < 1/ (12iVlog (NQ)). For given < r < 1/2 and < a < WNl f g(Ne) , the 
last inequality implies 

^h 2 e-W/ ia dx < N\og(Ne) [ h 2 dx + l§ [ \^-h 2 e^ 2 ' ia dx 
<Nlog(NQ) J h 2 dx + \J l^-h 2 e~^ 2 / 4a dx . 

Thus, 
(2.10) 

J l£h 2 e-W 2 / 4a dx < N log (NQ) J h 2 dx when 0<a< 16jvlo r g 2 (jve) ■ 

Now, e-W 2 / 4a \x\ 2 /(16a) > (NO)~ N N log (NQ) when r < |x| < 2r and a = 
i6Anog(jve) ■ This and 1 2 - 1Q P impiy 

/ /i 2 dx < (NQ) N I h 2 dx when < r < 1/2 . 

□ 

Proof of Theorems® and [3j Constant coefficients. Let u satisfy (|2.2|) and set / = 
uip in Lemma [21 where -0 £ C^fij), 0<^<1, ^> = lin-B3 and ^ = outside 
Br. Then, 

2 

(2.11) |A/ + d t f\ < N [|/| + |V/| + \\u\\ LHQi) x Bi \B 3 ] 

and from Lemmas [1] [2] and (|2.11[) , we have 

H a {t) > N-\t + a)- n ' 2 e- 1/(t+a) f u 2 (x, 0) dx , 



(2.12) JB i 

^a(t) > -n \i + N a (t) + e - 5/4{t+a) e 

when < t+a < l/Nlog(NQ) and where was defined in Lcmma[TJ The choice in 
Lemma[T]of the constant N gives, e - 5/4(t+a) e < 1 when < t + a < l/N\og(NQ). 
This fact and the second inequality in (|2. 12|) imply that 

(2.13) e Nt N a (t) + Ne Nt is non-decreasing when t + a< l/Nlog(NQ) . 

Multiplying (|2.5[) by H + (t) > we nnd f rom (EH]), Lemma [T] and the first inequality 
in (|2.12p that there is some constant N — N(n, M) such that 

(2.14) N a (t) < N [1 + (t + a)d t log H a (t)] , when < t + a < l/Nlog(NQ) 
and 

(2.15) d t lagH a (t) <N[l + N a (t)] /{t + a) , when 0<t+a< l/Nlog{NQ) . 

Setting [3 = jviog^wo) ' we f rom <|2-13[) . (|2. 14[) and the first inequality in 
(p2~T2f that 

r/S/2 N ( t ) 

N a (t) < N a {(3/4) + 1 < 1 + / dt<l+ d t logH a (t) dt 

(2.16) h/i ( t + a ) J p/4, 

= 1 + log < Nlog (NO) , when t + a < (3/12 
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and with constants depending only on M and n. In particular, the inequality 
2a J |V/(x,0)| 2 e~ |K|2/4a dx+f / f(x, 0) e -^ 2 / 4a dx 

<Nlog(N@) ( f 2 (x,0)e-W 2/4a dx 



(2.17) 



holds for all < a < 1/ (12AHog (NO)) 

From Lemma S] with h = /( -,0) and recalling that f(x, 0) = u(x, 0) in B3 we 
obtain the doubling property claimed in Theorem [5] when z — 

(2.18) / u 2 (x,0)dx < {NQ) N / w 2 (x, 0) when < r < 1/2 . 

For the same values of r choose k > 2 such that 2~ fc < r < 2~ k+1 and iterate 
([2~18| when r = 2~*, j = 0, . . . , k - 1. It gives 



M 2 ' 



! (a;,0)dx < (7V9) wlog(1/r) / M 2 (x,0)dx , 

but writing the value of in the above inequality, which was defined in Lemma [T] 
one gets 

1 JV log (l/Q 

/ /■ „ \ l + JVlo E (l/0 / f „ \ l + «log(l/0 

u 2 {x,0)dx < i^j u 2 (x,0)dxj i^N j u 2 dX j 

which proves the two-sphere one-cylinder inequality. 

Rewrite the integral on the left hand side of the following inequalities as the 
sum of the corresponding integrals over B r and R™ \ B r and assume that a < 
r 2 /(16AHog(iVe)), then 

f 2 (x,0)e-^ 2 / 4a dx< [ f 2 ( x ,0)dx + ^ [ \£.f{ x ^) e -^ 2 / 4a dx 

J B r JR"\B r 

< f B f 2 (x,0)dx + m ^ m J i£f(x,0)e-^dx 

and recalling that ([2"TU|) holds when h = /( -,0) and a < r 2 / (16iVlog (iV6)), it 
follows that 

(2.19) [ f 2 (x,0)e-W 2/ia dx < 2 [ u 2 (x,0)dx when < a < 



g(»e) ■ 

The fact that f(x,0) = u(x,0) when \x\ < 2 and choosing a = lgjv i^-pygj in 
the last inequality and (|2 . 1 7[) implies 



(2.20) / \Vu(x,0)\ 2 dx < (NS) N r- 2 / u 2 {x,0)dx 



|2 j„ ^ fATdk\N„-2 I „,2, 

Then, (|2.20| and the Sobolev inequality 



give 



Tl-2 1 

-j- |u(x,0)|^ dx\ " <{NQ) N j- u 2 (x,0)dx , 
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which shows that |u| 2 |f = o satisfies a reverse Holder inequality centered at 0. When 
n = 2, replace the previous Sobolev inequality by 

. i i 
-j. \if-~j- ipfdxj < Nr (J- \V<fi\ 2 dxY , ipeC°°{B r ) . 

These and other standard arguments show that |it| 2 |* = o is a Muckenhoupt weight 
in B 1 / 2 , which satisfies a reverse holder inequality from L™- 2 -averages to L 1 - 
averages and with constant 



N 



\lB 1/2 u2 ( X '°) dx J 

To prove Theorem El (EUSJ) and ([2TT5]) imply 

d t \ogH a {t) < Nlog(NQ)/(t + a) , when < t + a < l/N\og(NQ) 
and integration of this inequality over [0, t] gives 

when < t + a < 12 jvi g(jve) • From (|2.19p and the last inequality 
(2.21) f u 2 (x,t)dx<(NQ) N (l + ±) Nl ° s{Ne) [ u 2 (x,0)dx, 

when < t < 4r 2 < 32JV i g (Jve) and a = 16Nl f g{Ney The integration of (|2.21|) 
over [0, 4r 2 ] gives 

\ N „ f „2^nu, r /Q n j_ ^iviog(ive) 



u z dX < (NO)" a / u'{x,0)dx / (1 + *)""*^ ds . 

Qir JB r JO 

Thus, 

(2.22) / u 2 dX <e NlosiNe ^°^ Nlo ^ N6 ^r 2 f u 2 (x,0)dx . 

jQ 2r JB r 

The second part in Theorem [3] is a consequence of (|2.22[) . the doubling property 
in Theorem [5] and the standard local L°°-bounds of solutions of (|2.2p in terms of 
their L 2 -averages [12]. □ 



3. Variable coefficients 

In this section P is a backward parabolic operator verifying (| 1 . 2[) and (|1.3[) and 
u a solution of the inequality 

(3.1) \Pu\ < M{\u\ + |Vu|) in Q 4 ■ 

As in section [5] we need a few Lemmas to prove Theorems [5] and |3] The first one, 
Lemma [5] is the variable coefficients version of Lemma [T] The second, Lemma HI 
is a more detailed version of the Carleman inequalities established in [9] and [ID] . 
These two imply Lemma [7] a result analogous to (|2.17p . which as we saw in section 
[2 implies Theorems [2] and [3] 
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Lemma 5. There is N — N(X,M,n) such that the inequalities 

N\og{NVL p )>l and N u 2 (x,t)dx> / u 2 (x,0)dx 



J E>2 P J Bp 

hold when < t < p 2 / (Nlog(Nfl p )) and < p < 1. Here, 

J Qi u 2 dX 
Qp = p 2 J Bp u 2 {x,0)dx ' 

Proof. When p = 1 set as in Lemma [1]/ = uip, where (f € Cg^-E^), < (p < 1 and 
ip = 1 in B 3 / 2 - Define 



H(t) = J f'(x,t)G(x,t;y,0)dx , 

where now G{x, t; y, s) is the fundamental solution in R n+1 of the parabolic operator 
P*, the adjoint of P, y £ B\ and P*G(x, t; y, s) = —S(y tS ) is the Dirac operator at 

The Gaussian bounds for the fundamental solution [1] 

N^ 1 ^- s)-%e~ Nl * I ° l < G(x,t;y,s) < N(t- s) _ ^e"W^ 



and the fact that (|2.ip is satisfied in this setting [12 , show that the arguments in 
the proof of Lemma [1] can be repeated again, which proves Lemma [5] when p = 1 . 
The case < p < 1 follows rescaling to the case p = 1 and from the observation 



u 2 dX < / u 2 dX 



□ 



In Lemma [6] we use the following notation: for a given function a defined on 
some interval and a > 0, a a (t) = cr(t + a), G{x,t) — £~™/ 2 e -M / 4t j s the Gauss 
kernel and G a {x,t) = G(x,t + a). 

Lemma 6. Assume g v (0,0) — Sij. Then, there is a constant N = N(X,M,n) 
verifying the following property: for each number a > 2 there is a nondecreasing 
function a : (0, +oo) — ► M + satisfying t/N < a(t) < t when < at < 4 and such 
that the inequality 

a 2 J a~ a v 2 G a dX + a J a\- a \Vv\ 2 G a dX 

<N [ a l - a \Pv\ 2 G a dX + N ot a a sup / v 2 + \\7v\ 2 dx 
J t>o J 

(a/N) I \\7v(x,0)\ 2 G(x,a)dx + Na [ v 2 {x, 0)G(x, a) dx 



+a(a)- a 

holds for all0<a< ± and v e C5 (R™ x [0, f )). 

The proof of Lemma [S] is sketched at the end of this section. We now prove 
LemmaH where as in (|2Tf|) . / = uip with ip e Cg°(B i ), < ip < 1, i/j = 1 in B 3 
and ip — outside Bt_ . 
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Lemma 7. There are constants N = N(X,M,n) > 1 and p = p(X,M,n) in (0,1) 
such that 

2a J \Vf(x,0)\ 2 e~^ 2/4a dx+% J f 2 (x, 0)e~ |x|2/4a alx 

<N\og{NQp) J f 2 (x,0)e-^^ ia dx , 
when < a < 1/(12^ log (NQ p )) and N log (NO p ) > 1. Here, 



dX/ / u 2 (x,0)dx . 



Proof. For fixed a > 2 and a < — take as v in the Carleman inequality the function 
v = f(x, t)(p(t) — u(x, t)ip(x)ip(t), where ip = 1 when t < — and cp = when t > — . 
From (13TTT) and JO) we have 



|«| + |V«| + ||u||z,2(Q 4 )XB 4 x[0,f]\B 3 x[0 i] 



iPul < N 



the observation that <j\ a G a < JV Q+ 2a a+ ! 1 in the region B 



4 X 



(3.2) 
Then, 

[0, -] \ B3 x [0, — ], standard arguments used with Carleman inequalities and the 
fact that t/N < a(t) < t on (0, -), imply that the inequality 

(t + a)- a u 2 e -|x| 2 /4(t+a) dx < N a a a \\u\\ 2 L ^ Qi) 



(3.3) 



J B 2 



(a/N) / |V/(a;,0)| 2 e- |x|2/4a dx + iVa / f 2 {x, 0) e - |a:|2/4a dx 



holds when a > § + 2 and < a < ~ . 

For < p < 1, which will be chosen later, let N be a fixed constant larger than 
the ones appearing in Lemma [5] and ()3.3j) and assume that a > N log (NQp) and 

< a < Then, 
(3.4) 



(t + ay a u 2 e-^ 2 / 4 ^ dX>a 2 



JB 2 



> 



N 



ds 



u 2 (x,0)dx > 2? 



> 



J B 



1 



(t + a)~ a e~^u 2 dX 
u 2 (x, 0) dx 



ds 



2N \pe 



2a 



i 2 (Q4) 



n, 



The inequalities (|3. 3|) and (|3.4p show that to make sure that the left hand side 
of (|3.3[) is greater than four times the first term on right hand side of (|3.3p when 

2 

a > AHog (NQ p ) and < a < suffices to know that 



(3.5) 



1 



> 



8iV (iVe 2 )"^ 



Choose then p as the solution of the equation = 8Ne 2 . Then, ()3.5() holds 
when 8 Q_1 > NQ p , Thus, there are fixed constants, p = p(X, M, n) in (0, 1) and 
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IB, 



N = N(X, M, n)>l, such that 

(t + a)~ a u 2 e -l^l 2 /4( t+a) dx + N a a a \\uf L2(Qi) 

-(a/N) J |V/(x,0)| 2 e-^ |2/4a dx + Na J f 2 (x, Q) e -^ /ia dx 

2 

when a > iVlog (ND, p ) and < a < . In particular, the inequality 
2a f \Vf(x,0)\ 2 e~^ 2/4a dx + % J f 2 {x, Q) e -^ 2 / 4a dx 

<Nlog{NQ p ) [ f 2 (x,0)e-^ 2/ia dx 



(3.6) 



is true when < a < p 2 / (12AHog (NQ p )) and N\og(Nn p ) > 1. Recall that 
fl p = P / p 2 and rename the fraction N / p 2 as N in the above claim. Then, Lemma 
follows. □ 

Proof of Theorems® and [3j Variable coefficients. Lemmas |4] and [7] imply that 

(3.7) f u 2 (x,0)dx < {NO p ) N [ u 2 (x,0)dx when < r < 1/2 
JB 2r J Br- 
and (|3.7p and the argument used in section [2] to prove the interpolation inequality 
or two-sphere one-cylinder inequality (2) in Theorem[2j show that the interpolation 
inequality 

f / f \ l + «log(p/r) / f \ l + Nlog(p/r) 

(3.8) J u 2 {x 7 Q)dx < l^J u 2 (x,0)dxj J u 2 dX J 

holds when r < % . On the other hand, translations of u in the space variable 
show that the same interpolation inequality is true when we replace in (|3.8p the 
balls B p and B r by B p (x) and B r (x) respectively, and where x is any point in B\. 
Then, standard covering arguments combined with these interpolation inequalities 
at large scales f (|3.8[) and its translated analogues when r = p/2), show that we can 
find [3 in (0, 1) such that 

( ^ 1-/3 

2/™ n^ I / «.2f™ r\\ A„ \ at I „,2 . 



u 2 (x,0)dx<\ u 2 (x,0)dx\ [N u 2 dX 

Bi \Jb p J V JQt 

inequality which can be rewritten as (NQ p ) < NQ. This fact, (|3~7| and the 
arguments in section [5] finish the prove of Theorem 

Following the notation in Lemma [H the arguments in section [2] give that the 
claims 3 and 4 in Theorem [3] follow if we know that 

(3.9) d f \ogH a {t) < Nlog{N®)/(t + a) when < t + a < l/N\og(N&) . 

On the other hand, translations of u in the time variable, Lemma[5]and Lemma[7](In 
particular, applying to u(x,t+s) the arguments leading to the proof of (|3.6p . which 
satisfies (|3.1[) for some backward parabolic operator P when s < p 2 /N\og(Nil p ), 
then recalling that in this case f(x,0) = u(x, s)tp(x) and then, rewriting the vari- 
ables a and saso + t and t respectively) , show that with a possibly a larger N and 
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a smaller p 

(a + t) f \Vf{x,t)\ 2 e~^ /i{a+t) dx<N\og(Ntt p ) \ f 2 {x,t)e^ /4{a+t) dx 



when < t + a < p 2 /Nlog(Nil p ). The last inequality, Lemma[3]and the inequality 
(NSpf < NO imply that 

N a (t) < N log (NO) , 
J ^- t f 2 (x,t)G a dx<N\og(NQ) J f 2 (x,t)G a dx , 

when < t + a < l/Alog(A6). 

Only using the Lipschitz regularity in the space- variable of the matrix of coeffi- 
cients of P to bound the first order term arising in the formal calculation of P*G a , 
gives that \P*G a \ < N(\x\ 2 + t + a)/(t + a) 2 . This fact, (|3~TU|) and a calculation sim- 
ilar to the one leading to (|2.5[) , which was done to compute H a {t) in terms of D a (t), 
but which on this occasion is done to bound H a (t) in terms of D a (t) but replacing 
the backward heat operator by P, imply that (|3. 9[) holds and prove Theorem [H □ 

Proof of Lemma [BJ Here we do not give a complete proof of Lemma |H] since it is 
basically already done in [5] and [TO]- We only outline the main details. 

As in [9], if P is a backward parabolic operator satisfying \1.2\ , (|1.3|l and 
<7 y (0,0) = Sij, we consider the time independent backward parabolic operator Q 



d i (g ij {x)d j ) + d t where g ij (x) = g i3 (x, 0) , x £ 



In [9l Theorem 4] it is shown with an integration by parts argument that there 
are N = N(X, M, n) and So = 5q(X, M, n) such that the inequality 



(3.11) 



S- 2 J a 2 - a \D 2 v\ 2 GdX + a J a~ a ^v 2 G dX + J a i-"°hH \Vv\ 2 GdX 



is satisfied when v € C§°(M. n x (0, -)), 7 = a/5 2 , a > 2, < S < S and where a is 



<N J o-^lQvfG dX+N a j a+N J (v 2 +t\Vv\ 2 ) dX 
§)),7 = «/<5 2 , 

the solution of the ordinary differential equation 

ft N (&)] = ^r 1 . CT (°) - . *(o) = 1 . whcrc = ^ ( lo s f ) 1 

In fact, a(t) = /3(7t)/7 with 



/3(t) = £ exp 



The main point here is that for some TV, t/AT < (3(t) < t, when < t < 4 and 

(3.12) t/AT < cr(t) < t , when < t < 4/7 . 

If for fixed < a < — one repeats the calculations which led to p. lip in [S], 
but now working with v £ C yD (M. n x [0, -)) and carrying out the integration over 
R n x [a, +00), one must take into account the boundary terms which occur when 
integrating by parts with respect to the time-variable and add up on the on the 
right hand side of p. lip certain new terms. In fact, an analysis of the proof of 
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Theorem 4 in [9] shows that there is N = N(X, M, n), such that the sum of those 
boundary terms is bounded from above by 



(3.13) (r{a)- a 



(a/N) / \Vv(x,a)\ 2 G(x,a) dx + Na / v 2 (x, a)G(x, a) dx 



+N I l -^-v 2 (x,a)G(x,a)dx 



([3~T2f implies, ^- G(x, a)(j(a)~ a < N a j a+N ', when |a;| > 8, < a < 1. This fact, 
Lemma [3] and the writing of the third integral in (|3.13[) as the sum of the same 
integral over Bs and R™ \ Bs show that (|3.13p can be bounded by 



(3.14) a(a\ 



-(a/N) / \Vv(x,a)\ 2 G(x,a)dx + Not / v 2 (x, a)G(x, a) dx 



+N a 1 



a a+N 



v 2 (x, a) dx , 



when 5 is sufficiently small. Thus, there is N = N(X, M, n) and 6q = 6(M, A, n) 
such that if v G C^°(M™ x [0, |)), a > 2, < a < i and S < S , then 



,2-ai n2„,|2 



-a9(7t)„,2 



v z GdX 



1-aBM) 



Vv\ 2 GdX 



<N I a 1 - a \Qv\ 2 GdX + N a 1 a+N sup [ v 2 + t\\7v\ 2 dx 

J t>a J 



-a(a)~ a 



-(a/N) / \Vv(x,a)\ z G(x,a)dx + Na v 2 (x,a)G(x,a) dx 



and where the integration in the integrals with respect to Lebesgue measure dX is 
done over R n x [a, +oo). 

If we plug in the last inequality the function v(x,t — a), consider the change of 
variables t = s + a and observe that > j when < jt < 4, we get upon 

renaming the new variable s as i that the inequality 

<T 2 J a 2 a - a \D 2 v\ 2 G a dX + a 2 J <j- a v 2 G a dX + a J al~ a \Vv\ 2 G a dX 
<N I a^ Q |Q(w)| 2 G Q dX + ^ Q 7 Q+A, sup / v 2 + (t + a)\Vv\ 2 dx 

+a(a)- a 

is satisfied when a > 2, v G C^°(M™ x [0, |)), < a < ^, S < S , §o is sufficiently 
small, and where the integration with respect to Lebesgue measure dX is done over 
R"x[0,+oo). 

On the other hand, \g % i(x,t) — g l ^(x)\ < My/i. Then, choosing 8 sufficiently 
small it is possible to replace Q by P on the right hand side of the above inequality, 
which finishes the proof of Lemma [6] □ 



(a/N) / \Vv(x,0)\ 2 G(x,a)dx + Na / v 2 (x, 0)G(x, a) dx 



DOUBLING PROPERTIES OF CALORIC FUNCTIONS 



15 



References 

[1] G. Alcssandrini, S. Vessella, Remark on the strong unique continuation property for parabolic 

equations, Proc. Amcr. Math. Soc. 132 n.2 (2003), 499-501. 
[2] G. Alessandrini, E. Beretta, E. Rosset, S. Vessella, Optimal stability for inverse elliptic 

boundary value problems with unknown boundary, Ann. Scuola Norm. Sup. Pisa, CI. Sci. (4) 

XXIX n.2 (2001), 755-806. 
[3] G. Alcssandrini, A. Morassi, E. Rosset, Size estimates, Inverse Problems: Theory and Ap- 
plications, Contemporary Mathematics (G. Alcssandrini, G. Uhlmann eds.), vol. 333, Amcr. 

Math. Soc. 2003, pp. 1-34. 
[4] D.G. Aronson, Bounds for the fundamental solution of parabolic equations, Bull, of the Amer. 

Math. Soc. 73 (1967), 890-896. 
[5] B. Canuto, E. Rosset, S. Vessella Quantitative estimates of unique continuation for parabolic 

equations and inverse initial- boundary value problems with unknown boundaries, Trans. Am. 

Math. Soc. 354 n.2 (2002), 491-535. 
[6] X.Y. Chen, A strong unique continuation theorem for parabolic equations, Math. Ann. 311 

(1996), 603-630. 

[7] L. Escauriaza, Carleman inequalities and the heat operator, Duke Math. J. 104, n.l (2000), 
113-127. 

[8] L. Escauriaza, L. Vega, Carleman inequalities and the heat operator II, Indiana U. Math. J. 

50, n.3 (2001), 1149-1169. 
[9] L. Escauriaza, F.J. Fernandez, Unique continuation for parabolic operators, Ark. Mat. 41 

(2003), 35-60. 

[10] F.J. Fernandez, Unique continuation for parabolic operators II, Comm. Part. Din". Equat. 28 

n. 9 & 10 (2003), 1597-1604. 
[11] N. Garofalo and F.H. Lin, Monotonicity properties of variational integrals, A p -weights and 

unique continuation, Indiana Univ. Math. J. 35 (1986), 245—267. 
[12] O.A. Ladyzhcnskaya, V.A. Solonnikov, N. N. Uralceva, Linear and quasilinear equations of 

parabolic type, Translated by Smith, S.. American Mathematical Society. Providence, Rhode 

Island. 1968 

[13] E.M. Landis, O.A. Oleinik, Generalized analyticity and some related properties of solutions 
of elliptic and parabolic equations, Russian Math. Surv. 29 (1974), 195-212. 

[14] F.H. Lin, A uniqueness theorem for parabolic equations, Comm. Pure Appl. Math. 42 (1988), 
125-136. 

[15] K. Miller, Non-unique continuation for certain ode's in Hilbert space and for uniformly 
parabolic and elliptic equations in self-adjoint divergence form, Arch. Rat. Mcch. Anal. 54 
(1963), 105-117. 

[16] C.C. Poon, Unique continuation for parabolic equations, Comm. Part. Diff. Equat. Appl. 

Math. 21 (1996), 521-539. 
[17] C. Pucci, Alcune limitazioni per le soluzioni di equazioni paraboliche, Ann. Mat. Pura Appl. 

IV Scr. 48 (1959), 161-172. 
[18] J.C. Saut, E.Scheurer, Unique continuation for evolution equations, J. Differential Equations 

66 (1987), 118-137. 

[19] CD. Sogge, A unique continuation theorem for second order parabolic differential operators, 
Ark. Mat. 28 (1990), 159-182. 

(L. Escauriaza) Universidad del Pais Vasco / Euskal Herriko Unibertsitatea, Dpto. 
Matematicas, Apto. 644, 48080 Bilbao, Spain. 
E-mail address: mtpeszulSlg.ehu.es 

(F.J. Fernandez) Ceremade, Universite de Paris IX- Dauphine, UMR CNRS 7534, France. 
E-mail address: fernandezaceremade.dauphine.fr 



(S. Vessella) DiMaD, Universita di Firenze, Via Lombroso 6/17, 50134 Firenze, Italy. 
E-mail address: serglo.vessellaadmd.unifi.lt 



